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ABSTRACT 


\ 


A kinetic theory for fluid turbulence is developed from the Liouvllle 
equation and the associated BBGKY hierarchy. Real and Imaginary parts of 
Fourier coefficients of fluid variables play the roles of "particles." 
Closure is achieved by the assumption of negligible five-coefficient 
correlation functions. Probability distributions of Fourier coefficients * 
rather than moments* are the basic variables of the theory. Neglect of 
the correlation functions does not necessarily imply a cumulant discard 
hypothesis or a quasi-normal assumption* However* a later additional 
approximation leads to a closed moment description similar to the so-called 
eddy-damped Markovian approximation. A kinetic equation is derived for 
which conservation laws and an H-theorem can be rigorously established. In 
the absence of viscosity and external driving forces* the H-theorem Implies 
relaxation to the absolute equilibrium of Kraichnan. The equation can be 
cast in the Fokker-Planck form, and relaxation times estimated from its 
friction and diffusion coefficients. An undetermined parameter in the 
theory, about which significant assumptions have to be made* is the free 
decay time for triplet correlations. Some attention is given to the 
inclusion of viscous damping and external driving forces. 
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I. INTRODUCTION 




Two major and distinct areas of Investigation in non-equilibrium 
statistical mechanics over the last thirty years have been the BBGKY 
kinetic theory of gases 1 ’’* and the statistical theory of fluid turbulence 8 ” 11 . 
It seems worthwhile to attempt to unify the two by considering the latter 
from the point of view of the former. 

The purpose of this article is to present a framework for a systematic 
kinetic theory of turbulence originating from the Liouville equation for 
the Fourier coefficients of the fluid variables. The real and imaginary 
parts of these Fourier coefficients play the role, in a somewhat abstract 
way, that particle coordinates (positions and momenta) play in the BBGKY 
theory* Early formulations in this direction are due to Edwards and 
Herring 1 ** 15 , following general indications of Hopf . 15 

l ? 1 IQ 

A related approach of some importance is due to Lundgren * (see also 
19 20 

Monin and Chung ). Lundgren' 8 basic variables are configuration-space 

probabilities instead of Fourier-space distributions* Some apparent success 

21 22 21 * 

using Lundgren* s formulation is due to Fox, * * who has reproduced the 

decay of an energy spectrum for grid turbulence. The essential approximation 
in Fox's theory has been the neglect of the three point spatial correlation 
function, by analogy with the kinetic theory of plasmas. However, in the 
fluid context, the physical content of this approximation is quite obscure. 

The approximation is justified for the plasma ca3e* by arguments based on 
thermal equilibrium theory which are lacking in the case of fluids. 
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One argument for the use of Fourier coefficients as dynamical variables 

le that* at least In the absence of viscosity, there is an absolute equili- 

25 26 

brlum distribution both In two and three dimensions (see Kraichnan * 

12 

and Lee ). It can be said that relaxation to this absolute equilibrium 

can be considered to have been definitively established in two dimensions, 

27 28 

by the numerical investigations of Seyler et al, * as the state to 
which the truncated lnviscld Navier-Stokes system relaxes for arbitrary 
initial conditions. This equilibrium theory can be used as a guide for 
suggesting approximations for the non-equilibrium theory, as in the case 
of the BBGKY kinetic theory. No such equilibrium theory is available in 
the configuration-space representation. 

It is our goal to introduce only quantities that can be given a sharp 
definition in terms of Liouville's equation and probability distributions 
(ensembles) which obey it. It is clear that at one level, everything about 
the statistics of the system must be obtainable from Lio’- '.le’s equation. 

It is also our intent only to make approximations which appear to be in 
some well-defined sense experimentally testable; basically two will be 
required. 

The approach is illustrated by considering what is algebraically the 
simplest case: two-dimensional, inviscid, incompressible Navier-Stokes 
flow. The methods are not restricted to this case and can also be 
applied to three-dimensional flow, to magnetohydrodynamic turbulence, and 
to cases when dissipation and external driving forces are present. Our 
intent is to introduce the method as simply as possible; detailed compa- 
risons with data will have to await numerical solution of the kinetic 
equation to be derived, and additional refinements and revisions by other 
workers are to be expected. 
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The basic relations of the problem are the incompressible Navler-Stokes 


equations In two dimensions with zero viscosity: 

du 

H + « * 7 H " " ?<P/P M > » 

V ♦ u ■ 0, 

where u Is the fluid velocity , p is the pressure , and p M is the (uniform) 
density, u lies in a plane normal to the z-axis, and is a function of x and 
y only. It is convenient to work in the vorticity representation, for which 

A 

p ■ pb * Vxu. 


A 

b denotes a unit vector in the z-direction. All variables are represented 

as Fourier series in a large but finite box. Both u and p are eliminated 

in favor of an equation Involving p alone in the usual way (see Eq. (1) 

below). A well-known alternative interpretation of the equations, that 

29 

they describe the electrostatic guiding-center plasma, is possible, in 
which p is the electrostatic charge density and u is the E x B drift 
velocity of a fluid element. 

In Sec. II, a notation is Introduced which expresses neatly the 
dynamics of the real and imaginary parts of the Fourier coefficients in 
terms of a simple set of nonlinear coupled ordinary differential equations. 
A Llouville equation and BBGKY hierarchy can be derived in the phase 
space defined by these coefficients (Sec. III). Though the moments of 


5 


the probability distribution functions of the Fourier coefficients 
have been more popular as a set of dynamical variables for statistical 
theories of turbulence, our concern throughout will be with the distri- 
bution functions themselves. 

Closure of the hierarchy is achieved in Sec. IV by the hypothesis that 
the five— coefficient correlation function is negligible. This closure is 
explicitly shown to include all cumulants in the moment hierarchy, and 

30 

iKt to be a variant of the cumulant-discard , or quasi-normal, hypothesis. 

The theory at this point owes a debt to the direct-interaction approximation 
of Kraichnan, where the assumption is made that the three-coefficient 

correlations are more basic to the dynamics than any higher grouping. But 
we have not introduced the unit infinitesimal response matrix of Kraichnan, 
an d have limited our apparatus to quantities which can be defined in terms 
of the single-time Liouville distribution. The kinetic equation derived 
in Sec. V shares, however, the non-Markovian property of the direct-inter- 
action approximation: The time derivatives of the distribution functions 

depend not only upon their present values, but on their recent history 
over a time of the order of a decay time for triplet correlations. A 
second hypothesis (for which experimental evidence for or against appears 
to be lacking) will result in a Markovian kinetic equation: That the decay 

of three-coefficient, correlations in the absence of their source terms 
is rapid compared to the evolution of the one-coefficient distributions. 
Numerous features of this Markovian kinetic equation (Sec. VI) can be 
rigorously established, including that of a monotonlc approach to the 
absolute equilibrium distribution of Kraichnan • (H-theorem) , a result 
believed to be the first of its kind. With no further approximation. 
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the kinetic equation can be shown to be of a standard* (nonlinear) 
Fokker-Planck type. The Fokker-Planck equation is of the kind derived 

OR 

from the Langevln equation for Brownian motion, with the difference 
that the friction and diffusion coefficients are time-dependent Integrals 
of the distributions of the other Fourier coefficients. At this level, a 
closed description in terms of low-order moments is possible, but the 
Fokker-Planck equation provides information on moments of arbitrarily 
high order. A relaxation time can be estimated from the Fokker-Planck 
equation, to be compared with the decay times for triplet correlations, 
when the latter become available. 

Only two significant hypotheses or approximations (beyond the standard 
but significant one of the representation of the fluid variables as truncated 
Fourier series) are employed, and both involve only well-defined quantities 
for the system: 

(1) Neglect of the quintuplet correlation and a subsequent assumption 
concerning the mode of relaxation of the triplet correlation; and (2) the 
assumption that the decay of the triplet correlation is rapid compared to 
the evoiw on of the single-coefficient dist ibutions. A significant 
fraction of the theory (specifically, everything up to Eq. (32)) can be 
based on the first hypothesis alone, and we re-emphasize that thiu is not 
quite a cumulant discard hypothesis. The second assumption is analogous 
to the Bogolyubov^ "functional assumption." Both are subject to some 
degree of experimental verification or disproof. 

Section VII discusses the problems associated with the inclusion of 
viscosity and external driving forces. Section VIII summar 4 ss the results 
and anticipates future inquiries. 
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II. BASIC DYNAMICAL DESCRIPTION 


For our purposes, it is useful to write the Fourier- transformed, 
two-dimensional, inviscid. Navier-Stokes equation in the vorticity 
representation: 


3p (k.t) m 

M(r,g) p (r.t) p (g.t) , 

p+r»k 

W ^ 


( 1 ) 


where the coupling coefficients between the various Fourier modes are given 
by 


b • (p x r) . . 

M(r,p) = M(p,r) « f- “o 

-v ~ t r p 



If the fluid velocity is u » u(x,t), the vorticity vector is p(x,t) » 

^ ^ ^ 

A A 

Vxu » bp(x.t) . where b is a unit vector lying in the z-direction, say. 
p(x,t) is written as a Fourier series 


p(x,t) - T. p(k.t) exp (ik • x) 

and p(k.t) « L -2 fix p(x,t) exp(-ik • x), and the integral runs over a large 

square box of edge L. The vectors u, x, k, r, p oil lie in the xy plane. 

' ***««**#*» 

A 

and are normal to b. The allowed values of the wave numbers k. p, r are 

2irn/L» where n is a non-null vector (n ,n ) with integer components. The 

** ^ x y 

wave numbers' magnitudes are bounded from below by virtue of the finite size 
of the box and from above by a large but finite maximum, which remains 
somewhat arbitrary. Extensive numerical studies of the relaxation of Eqs. (1) 
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were carried out by Seyler et al, starting from highly non-thermal initial 
conditions . 

It Is the statistical theory of Eqs. (1) which is of interest, but 
it is convenient to Introduce a more compact notation (see, e.g. , 
Kraichnan,^ Betchov,*^ or Herring^). We order the wave numbers by 
associating them in any convenient way with the integers and represent them 
by Roman subscripts i» j» k, 1, ... . There are to be two integers per 

wave number. By the variables X^, we shall mean the real 

and imaginary parts of the various p (k,t) associated with the ordered wave 
numbers. 

The system of equations represented by Eq. (1) can then be written in 
the simple-looking form 

d x i v 

-3 r - £ c 1Jk X J ^ . c 


where the are a set of real constant coefficients which can be Inferred 
in terms of the M(r,p) from Eq. (1). All terms in Eq. (3) are now real. 

The subscripts range from 1 to N, say, where N»l. N is the number of 
allowed wave numbers. (There are 2N total real and imaginary parts of 
Fourier coefficients, but only half the variables are independent, since 
p(k,t) » p*(-k,t).) Only variables which can be prescribed: independently 
initially are included in Eq. (3). 

For detailed numerical applications, the exact values of the co- 

efficients can be read off, but for the formal theory to be presented here, 
only five properties of the are needed: 


(X) .'ll c 


0 unless 1, 4» k correspond to different wave 


ljk " 

numbers (l.e. , unless l f k are all different)! 


(2> C ljk" C ikj { 

(3) 1, J, end k all run from 1 to N; 

M) C 1Jk + C JU + c ktj - 0; 

(■» £iis_ + fjM + 5sii . 0 

p J 1 ' 2 Ul 2 w 2 

Property (A) expresses conservation of enstrophy (or mean-square charge 
density, for the plasma case). Property (5) expresses conservation of 
energy. In property (5),|i| 2 Is a symbolic notation which means the square 
of the wave number associated with the subscript 1. (The constants of the 
motion guaranteed by properties (4) and (5) survive an essentially arbitrary 
Fourier truncation.) 


Equations (3) may also be used to discuss the three-dimensional case, 
and much of the theory to be developed applies to the three-dimensional 
case as well. The major difference is that in three dimensions, the 
analogue of property (4) does not apply. The consequences that result if 
property (4) does not hold will be noted later. 

The initial-value problem involves specifying the initial values of 
the Xj and following their subsequent evolution according to Eqs. (3). 

We turn now to the introduction of probabilistically distributed initial 
conditions. 


III. LIOUVILLE'S EQUATION AND THE BBGKY HIERARCHY 


The various In Eq. (3) con be thought of a coordinates In a 

many-dimensional phase space. Eq. (3) apparently does not define a 

Hamiltonian system. It Is* however, a conservative system and possesses 
25 

two constants of the motion as a consequence of properties (4) and (5) 
of Section II: The energy £ and enstrophy ft; these are 

£ - E t X*/|l| 2 (4) 

“ ■ (5) 

Moreover, a Liouvllle equation, or probability conservation law, exists and 

governs the evolution of any probability distribution law D - D(X. X„ ... X ,t) 

x 2 n 

defined over the phase space: 

«♦ ; •) - °- «> 

In Eq. (6), dXj/dt is given by Eq. (3). Using property (1), 

t? + £ c i Jk *1 h 557 » - “ • 

Conservation of the positive semi-definite character of D and conservation 
of the normalization integral 

/ dX £ . . . dXj, D - 1 
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" r ■¥»'- T ir-pf pjr ^ ft ffif j »;' »v» fwrnw 7 jn |nw pq* hipri ww w * in irtwff*'' m 


'V Will upfipm WW wpiwrw-wTT 


rr ■ r- • • ^««iiiuiii p-nnim 


follow immediately from Eq. (7). 

A family of redueed probability distributions f... - f... (x. , 

*4 K* • • r lj R • • ♦ r 1 

Xj* ••• X r ,t) can be defined by Integrating D over all the phene space 
coordinates which are not of Interest: 


s /»(« 1«2 •” *h) 


ijk. . ,r 


(8) 


where the notation in the integral means to integrate D over all the phase 

space variables except those associated with the indices ijk...r. The reduced 

distribution functions defined by Eq. (8) differ from those used in the 

1-7 

conventional BBGKY theory mainly in that they are not symmetric under 
Interchanges of the different phase space coordinates. Thus f ^(X^,t) is not 
the same function of X i that fj(Xj,t) is of X jt if i and j are different. Most 
of the added complexity of the theory sterna from this lack of symmetry. Note 
that the subscripts on the fjjj^ r are by definition always different. There 
are never any repeated subscripts in the distributions. 

The basic variables in most statistical theories of turbulence are 
moments, such as 


<lt i > 5 I A.* i dx i 


or 

<x i y s fh hi dX i dx j • 


However, in these present calculations, the variables to be treated as basic 
are the distributions f^, f ^ , f^ etc. 
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By integrating Eq. (7), discarding surface integrals at infinity in 
phase space, and using Eq. <8), a BBGKY hierarchy of equations for the 
reduced distributions can be readily generated. For example, the first two 

members are 


and 


" "db L C ijk / dX J dX k X J *k f ijk 
1 jk 


(9) 


-5? + ( c tjk x j *k + '• »• ) f i. k 

’ - IT £ C lr. / dX . dV ' . ijkrs + - »• 
i rs 


( 10 ) 


The notation "c. p." will always mean, add tc the immediately preceding term 
the samf. term with the unsummed indices cyclically permuted, pud continue 
doing this until the original ordering is reached. Thus 


C ijk X j \ 3X t + C * p * 

C ijk X j \ 3X^ + C jki C kij X i X j * 

As in the case of the BBGKY hierarchy for interacting point particles, 
the higher members of the hierarchy become complicated. So far, there 
has been no incentive to define a notation complicated enough to write 
down the general term. But its structure is obvious. For s phase space 
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coordinates, it will Involve the Liouvllle equation for s coordinates on 
the left and a large number of similar integro-dif f erentlal linear 
operators on the right, which act on distributions of s+2 coordinates. 

One interesting feature is that reduced distribution functions containing 
only odd numbers of the X^'e directly enter the hierarchy. Though the 
distributions of even numbers of coordinates can be obtained from the 
odd distributions by integration, they play no direct dynamical role. This 
« ontrasts with the all-important role that the pair correlation function 
plays in particle kinetic theory, and with Lundgren’s configuration-space 

17 

description. 

A cluster expansion can be used to replace the multibody distributions 
f ijk* £ ijkrs* • * * » b y appropriately expressed correlation functions 
Qijkrs* •••» ®® follows: 

( 12 ) 


( 13 ) 

The "quintuplet" correlation Q^j^rs *® tbe completely nonfactorable part 
of the and so forth. 

We shall restrict attention to the case of vanishing first moments 

<X i > 5 / X ' f i dx t " °» all i. ( 1 A) 


'ijk " f i f j £ k + T ijk 


ijkrs 


f i f l f k f r f . + f l { i T kr. + f l f k V 

+ f l f r V + f l f . T jkr * f J f k T lrs 

+ { i f r T iks + { i *s T ikr + f k f r T «. 

+ f k f . T ijr + f r f s T ljk + ^ijkrs - 
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This corresponds to the conditions <p> * 0 t <u> ■ Q* etc# in conflgu- 

ration space.® However, we do not impose all the conditions that have been 

sometimes associated with more restrictive definitions of "homogeneous 

2 2 2 2 

turbulence." Also, it is not assumed that when i and 

j are identified with different wave numbers, and we do not assume that 
the initial values of the univariate distributions are Gaussian. 

From Eq. (14) , it follows that an alternative way of writing Eq. (9) 
is 


/ dx l dx k X j \ V • 

From Eq. (9) or Eq. (15), it is clear that what is required for a 
kinetic description, in the kinetic theory sense 5 * 6 ’ 7 of the term, is an 
approximate expression for the f or in terms of the f It is 

clear that some additional complexity (over and above that encountered in 
hierarchy derivations of Boltzmann's equation or the Fokker-Planck equation) 
is to be expected, since the different f^'s are distinct functions. 

So far, everything that has been said is exact. Further progress 
depends upon developing satisfactory approximations with which to close 
or truncate the hierarchy. 
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IV. CLOSURE OF THE HIERARCHY 


A hierarchy generated by taking moments of Eq. (1) or E~ , (3) has 
In the past been used as a starting point for statistical theories of 

Q 1A 

turbulence. * Closure schemes for this moment hierarchy, based on dis- 
card of higher cumulants,^ have largely fallen into disrepute. 9,10,37 

2 

This is in large part due to the negative spectral densities (<X^> , in our 
notation) which result. It is not true, as we shall later prove explicitly, 
that closure schemes for the just-derived BBGKY hierarchy based on the 
neglect of one of the higher-order correlation functions (such as 
are equivalent to the neglect of higher cumulants in the moment hierarchy. 

It will also be seen that neglect of the higher-order correlation functions 
need not, in general, lead to negative values of any intrinsically non- 
negative quantities such as <X^>. 

Host closure schemes for the interacting-particle BBGKY hierarchy 
rest on the identification of some small parameter which enters as a 
multiplicative factor in the hierarchy (e.g. , the density for Boltzmann's 
equation or the coupling constant for the Fokker-Planck equation) . This 
small parameter then serves as the basis for a perturbation expansion. ^ 

One of the more stubborn impediments in turbulence theory has been the 
lack of any such readily- Identifiable small parameters. The only para- 
meters which enter Eqs. (3), (7), (9), (10), or (15) are the coupling co- 
efficients C. and they can be considered in no sense "small." In 

xjK 

particular, the usual Fokker-Planck arguments (in which weak coupling is 
argued to imply weak correlations) would appear to be Inappropriate. 
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Nevertheless, there is a sense in which the correlations may still 
be argued to be small. One exact solution to Eq. (7) is the absolute 
equilibrium^* ^ distribution of Kraichnan: 


D » c exp { -a £ -fKi ) 
eq. 


where c is a normalizing constant, a 8 are constants which play the role 

of inverse temperatures, and £ and Q are given by Eqs. (A) and (5). 

Equation (16) implies Maxwell distributions for the individual X^s, 

«t> ■ c i“ p J-(|^ +s k|- ai 

and implies also that all correlations vanish identically. Extensive 
numerical evidence accumulated by Seyler et al shows that in fact a 
quite definite relaxation of the proper tits of the system to those pre- 
dicted by Eqs. (16) and (17) takes place. This relaxation to equilibrium 

seems to be beyond dispute. 

If we are in some sense not too far from the absolute equilibrium, 

predicted by Eqs. (16) and (17), then it is clear that the correlation 

functions must be in some sense small, so that, for example, fjf j f k >>T ijk* 

* f f f f »n . . and so forth. Dynamical arguments to support this 

TJ k t s x ijkrs* 

conjecture are given later in this section. 

The only assumption required for closure at a non-trivial level is 

that the quintuplet correlation Q 1Jkrs ls negligible compared to the product 

f f T, . This is somewhat in the spirit of the "direct interaction 
i j krs 

approximation,^*”*^ which assumes that correlations among three modes are 
more important than among more elaborate groupings. 
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Dropping terms containing Qjj^* expressing f Krs in terros of 
the f^^ and the f^, using Eq. (14), and performing some tedious but 
entirely straightforward algebra on Eq. (7) reduces It to the following 
form 


at + c ijk x j *k ax A + c * p * 


ijk 


I C irs IT / dX r dX . X r X . f i 

I1» /«. 1.4 1 l> 


krs 


rs^jk or kj 


+ f l *k f jr. + f J f k f ir. 


j + c. p. 


’ W >ijk ■ <»> 

V denotes a linear differentio- integral operator containing the f ^ whose 
action on the f^ rg Is defined by the right-hand side of Eq. (18). 

Equations (9) or (15) and (18) constitute the basic closed description 
that results from the neglect of quintuplet correlations. It is not possible 
to proceed farther analytically without more insight into the properties of 
the linear operator V than is now available, tinless we make plausible 
simplifying assumptions about the effect of the linear operator V. This 
is deferred to Section V. Numerical solucion to the coupled pair of 
equations (9) and (18) may uot be significantly more difficult than the 
solution of the Eulerian direct interaction equations, but this has not 
yet been attempted. 

To dose this section, two other remarks are in order concerning 
Eqa. (9) and (18). 
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First, it can be seen easily that the neglect of the quintuplet 

30 

correlation is in no sense a M cumulant discard" approximation. For, 
applying fdX ± X* to Eq. (9) gives 


d_ 

dt 



L 

jk 


c ijk <X i X J 


V 


( 19 ) 


while applying UX l dXj dX R X ± Xj X R to Eq. (18) gives 


<L < x 

dt X 


i X J V 


C ijk <X 


54 


> + c. p. 


( 20 ) 


Instead of (20), we would have in it ha cumulant discard approximation 


3t <x i x j V ' c «k <x ^ > <4 +c ' 


( 21 ) 


It is clear that at the level of the moment hierarchy, the neglect of the 
quintuplet correlations does not even provide closure. Any moment-develop 
ment equation will always involve higher moments because of the three- 
coefficient Liouville operator on the left-hand side of Eq. (18) . 

Finally, the content of the relation (13) for f^ can be translated 
into an equation for using Eqs. (12) and (15). The result can be 

written compactly as 


!_ 

dt 


+ 





( 22 ) 


where S 


ijk 


is a source 


term involving only one-body distributions: 




+ c . p, 


f i f j f k ‘ 


(23) 
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The notation "c. p." has its usual meaning, and H is a linear operator 
defined by the following relation: 


H T 


ijk = 


l iJk 


C ijk X J \ 3X t + c * p# 

1,. C ir. /« r dX . X r X . ST ", T kr. + f k T ,r.> + «• »• 


rs/jk 
or kj 


- 2 l c «k ‘i f k dx j 




(24) 


(H T^ k involves all the various triplet correlations, and not just T^ k itself.) 

If ve like, we can choose to regard Eqs. (15) and (22) as the basic dynamical 
equations. Again, the only approximations involved in them are the neglect 
of the quintuplet correlation functions. If we make the conventional and 
reasonable assumption of initially-vanishing correlation functions, it will 
be seen that acts as a source term which generates triplet correlations 

and involves only the f^. Note also that if f ^ is a Maxwellian, as in 
Eq. (17), S ljk 5 0. Thus for situations close to absolute equilibrium, the 
triplet-correlation-generating term in Eq. (22) is small even though the 
Cijk are not sma ^^* This provides a dynamical argument for the weakness 
of the correlations that is independent of the weakness of the interaction 
coefficients. 

It will be noted that the general structure of Eq. (22) is reminiscent of 
the Dupree form of the equation for the pair correlation encountered when 

7 

deriving the Balescu-Lenard equation. 
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V. APPROXIMATE SOLUTION FOR T^^} THE KINETIC EQUATION FOR TURBULENCE 

We coni ec ture that the effect of the last two collections of terms In 
Eq. (24) is a large number of essentially random Impulses which tend to 
destroy correlations and compete against the source which 

attempts to build them up. The nature of the decay of correlations pro- 
duced by the last two brackets of terms In Eq. (24) Is undoubtedly very 
complicated , but it may not be too poorly represented by a simple exponential 
decay. That is, we conjecture that Eq. (22) can be represented by the 
approximate form 

jit + + v ijfc J T ijk " S ljk * 

In Eq. (25), Lj(i,j,k) is the Liouville operator for three coefficients, 
and Is defined by 


I* 3 (l,J.k) = C 1Jk Xj X R 3X + c. p. 


V Is an Inverse decay time for three-body correlations. For many pur- 
poses it is not necessary to assume a specific form for it. Roughly speaking, 
plays the role that the Landau damping decrement plays in the derivation 
of the Balescu-Lenard equation. 

The formal solution of Eq. (25) for T^ k is, assuming zero initial 
correlations. 
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1|'U-'V..| * '■ . -1| ’ «' V; * | '/i.-, .* i*|„ 1 » ’ .st ij ( ; « ■ i, i • <1 ’1 s l 


T ljk " / dT exp j’ (t “ T)L 3 (1 *** k) I ex P |” (t ~ T > V ijk| s ijk* (26) 

The correctness of the solution (26) can be readily established by 

differentiation. 

The exponential operator exp{-(t - t) L 3 (i,J,k)} Is a familiar operator 5 * 6 * 7 
In BBGKY kinetic theory, and Is a member of the class often called "streaming 
operators." Their effect on an arbitrary given function Is to treat the 
function as an initial value for Liouville*s equation and convert It Into 
that solution of Liouville's equation which reduces to the given function 
Initially. More specifically, for any arbitrary function g(X i , Xj, X^) , 

expt-T L 3 (i,j,k)] g(X ±t X j , X k ) 

- g (X 1 (-T), XjC-T), X k (-T)). : (27) 

In Eq. (27), X^(x), Xj(x), X^(x) are the solutions of the ordinary differential 
equations 

d X. (x) 

dX “ C ijk X j (T) 

d X (x) 

“df~ " C Jk i V T > V T > 

d \(T) 

“di~ " C kij X i (T > X j (T) < 28 > 

which satisfy X^O) - X ±t X^(0) - Xy X R (0) • X k< The X £ (x) are Lagrangtan 

coordinates which are associated with the Eulerlan coordinates X^ by the 
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differencial, equations (28) and their accompanying Initial data (X^O) * Xj_» 
etc.). The explicit solutions to Iqs. (28) are net difficult to extract, 
but their explicit values are not required for present purposes. 

Wth these notations and a slight change in the variable of temporal 
Integration, Eq. (26) becomes 


ijk 


/• 

» I dT exp j-T v 1Jk | S 1Jk (X^fT), Xj(+t), X^-rr), *) 


(29) 


Equation (29) Is our approximate solution for T^ k » Its substitution In 
Eq. (15) gives: 


3 f ± (X^t) 
_ 


A £ C « k / “J dX k X J *k / dT exp( ' T 'V 

X [c iJk X j( .T) X^.T) + =• >•] 

f ± (X 1 (fT)t^ fj (Xj (ft) .t^f k (X k (4T)i) 


Equation (30) Is our basic kinetic equation for turbulence. It can be 
further simplified, and a number of its properties can be proved. This 
Is done in Sec. VI. (A formally exact restatement of Eq. (30), which 
does not assume the exponential relaxation times v^ k , is obtained by 
simply replacing exp(-x v^ k ) by the exponential of the operator whose 
eigenvalue V^ k Is conjectured to be; this formal generalization is at 
present of no computational value.) 


(30) 




23 


VI. PROPERTIES OF THE KINETIC EQUATION FOR TURBULENCE 


Three conservation laws can be proved directly from Eq. (30) without 
further manipulations. 

^ Conservation of probability follows from the obvious conclusion 

that 


It /f i dX i 


/dX 


a function of X, 


JL I 

* ** X i I which 0 as 


- 0 . 


^ 2 ) Conservation of ens trophy is proved by considering the expression 


a 

It 


- 

“ °ijk f dX i dX j dX k X i X j \ f dT ex P(" 


T V ijk> 


c ijk V* T) X k (+T) a^(«r + - c - p 


x (X 1 (K)J^) fj (Xj(tT)fi)' f k (X k (+T)ff) . 


(31) 


Permuting the dummy indices i, j, k in Eq. (31) and adding the three equi- 
valent expressions together, then dividing by 3, gives the same expression 
as Eq, (31) with the first C^ k on the right-hand side replaced by 

3 (C iJk * C Jki + C kiJ*» whlch b y Property (A) of Section II, is zero for 
all i, J, k, 

(3) Conservation of energy is proved in the same way as conservation 
ct enstrophy, with ± (C + C + C ) being replaced by + 

!usi + fkii\ J V ur 

Ul 2 Ml 2 / 
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A fourth obvious property 1st 

(A) Solution bjr i* Maxwellian . Since the right-hand side of Eq. (30) 
vanishes Identically when f i is given by Eq. (17), It is eleor that the 
absolute equilibrium f ^ Is a time-independent solution of Eq. (30) . 

The kinetic equation (30) ns it stands Is not Markovian. The time 
derivative of the f ^ depends not only upon at the present Instant, but 
on Its Integral over a time of the order of the the decay time for 

freely decaying triplet correlations. A Markovian limit can be obtained 
by making the following assumption, which is the analogue of the 
Bogolyubov 5 "adiabatic hypothesis." We may assume that the free decay of 
triplet correlations is sufficiently rapid that neither the f ^ nor the 
three-body orbit variables of Eqs. (28) change appreciably in a time of 
the order of v’^. This hypothesis either is or is not true, but experi- 
mental or numerical data that would shed light on it would appear to be 
scarce. If we make it, we may freeze the l— dependent values of the vari— 

o 

ables in the integrand of Eq. (30) at their au s s a a t values, and perform 
the x-integrations to get: 

TT * 5^ £ ^ \ - 

+ C jki *k X i atj’ + C kij X i X j 3xj[ 

X f ± (X A ,t) fj (Xj »t) f k (\.t) . (32) 

From the Markovian form (32), it is immediately possible to prove 
the positive semi-definiteness of f^ and that all approach a Maxwellian 
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a* C "* 00 (H-theorem). It can alao be ahown that a moment closure has 
now boon achieved. 

<5) Poaltive-dof lnltoneaa of f^ t Suppose that f^ (X^.O) la every- 
where > 0, all 1, and first becomes negative for some X. » X. at t ■ t 

l lo o 

as a consequence of Eq. (32). We will show that a contradiction Is 
implied. For at t • t Q , the stated conditions imply: 


<i) 


(ii) 

3 *1 <X io' ‘o’ 

3 ‘o 

(iii) 

8 f l < X io- 'o> 

3X lo 

(iv) 

32 f 1 «!.• *o> 

a v 2 


0 


< 0 


- 0 


> 0 


Evaluating the right-hand side of Eq. (32) under these conditions gives 


3 f i < X io’ 

5T““ 


*o> 


E 

Jk 


'ilk 


ijk 


/ dx ] dx k x j *i 


2 >2 h «lo' 


3 X’ 


io 


*.> 


<Y '„> f k ( Y 


*.> 


which is manifestly positive* contradicting (ii). Thus no Intrinsically 

2 

positive moment* such as can ever become negative as a consequence 

of Eq. (32). 
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(6) H- theorem . The irreversibility fan be demonstrated by considering 


a alight generalization of Boltzmann's H function. 
k £ /<•*! f t »■> f t - £/ln f t dX t 

- -£/«*«,«„ b t *,* 


x C_, X 


ijk x j \ ax A + c * p * I f i f j 


’ - I £ 7^- / dx i dlt , < i)t k 
J ijk v ijk j 1 J K 


c »k N \ iii 


+ c. p 


a x. 


f i f j *k - ° • 


since the two { }'s are identical and the rest of the integrand is > 0. The 
equality sign holds if and only if the { } ■ 0, or 


: ijk x j *k ax t + c- p ‘ 


f i f j f k ’ 0 • 


which is satisfied by Eq. (17). Since E^ /dX^ f ^ in f ^ is bounded from 
below, this completes the proof of the H-theorem. All initial distributions 
must approach Eq. (17). 

At the level of Eq. (32), in contrast to Eq. (30), a closure in terms 

2 

of moments has been achieved. Applying /dX^ (X^/2) to Eq. (32) yields, 
after some manipulations. 


( 33 ) 


<X^> 

d V 

dt 2 


%i; | c «k <x j 2> + c m <x i > 


+ c k« <*]> 


2 

and we have already proved that <X^> can never go negative, though a similar 

proof can be giveu directly from Eq. (33) . A somewhat similar set of 

38 39 

equations Is due to Leith and Orszag. Except for the slightly different 

form and interpretation of the v^, the content of Eq. (33) is that of 

38 39 

the eddy-damped Markovian model. * Relaxation to equilibrium for 
Eqs. (33) is guaranteed by our H-theorem. 

Equation (32) can be given a somewhat more simple looking form in 
terms of moments: 


8 *t 
it 


ax. 


< x i h £ i> + i 


(q 


i f i> 


(34) 


where 


X, 5 £ 

A it. 




^ V <=£ + c kij .<$■] 


and 


5 Z <X?> <X*> >0 . 

Jk v ijk 3 K 


(35) 


(36) 


Both and q t are independent of all the phase space coordinates and depend 
only upon the time Equation (34) has the classical form of the Fokker-Planck 
equation derived from the Langevln equation for the case of one-dimensional 
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Brownian motion with (time-varying) friction and diffusion coefficients, 
-A^ and q^. 

A relaxation time for the mode can be estimated by f )" 1 • -A ^ 1 


where 


'i : j ^ [ c jk i <> + Su • 


Using thermal equilibrium estimates for the expectation values of the moments, 
this is approximately 


X ± = - E <X 2 > <xj> . 

jk v ljk 1 * 

If typical <X 2 > are just <X 2 >, this gives a relaxation time of 


ft " 1 S -l " 1 S V 

1 1 C 2 <X 2 > 2 N 2 

where C 2 is a "typical" C 2 jfc , and v is a "typical" v^. 

The condition for validity of the Markovian assumption is that this 
relaxation time be long compared to the triplet decay time, or 


V » I 

C 2 <X 2 > N 2 V 


2 2 2 2 
V » C <X*> N . 
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The estimate for (which undoubtedly depends upon the the <X^>, 

and M) which would validate or Invalidate the Inequality (38) Is 
lacking. 
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VII. VISCOSITY AND EXTERNAL FORCES 


Viscous damping is relatively straightforward. One adds a term 
~ V i X i t0 tight-hand side of Eq. (3), where " |i| V, and v is 
the kinematic viscosity. The system is of course no longer conservative, 
and approach to a Maxwellian (17) is no longer expected. The new terms 
are reflected in new viscous damping terms in Eqs. (9). (10). (15), and 
(18). In particular, the operator H in Eq. (24) has added to it a term 
( V i X i + c * P* ^ T ijk on t *' e tight-hand side, and each of Eqs. (28) 

has a damping term like -v^ X^(x) on the right-hand side. The three-body 
orbits defined by Eqs. (28) are no longer conservative. Again, further 
progress requires conjectures about relaxations of triplet correlations 
in order to write down an explicit expression for a kinetic equation such 
..as Eq. (32). If we again assume that these triplet correlations' free 
relaxation rate is faster than any other characteristic time in the system, 
the net effect of the addition of viscosity is to add a term 3(v^ X^ f^)/3 X^ 
to the right-hand side of Eq. (32). The relaxation process then becomes 
a complicated competition between the non-dissipative "collision term" in 
Eq. (32) with the additional viscous term which leads to the decay of the 
energy and the concentration of the distributions at lower values of X^. 
Numerical investigation would be required. 

External forces are more difficult to add. The Inclusion of a driving 
force F^t) on the right-hand side of Eq. (3) will, if random, affect the 
evolution of D in a non-Markovian way. Simplifying assumptions which would 
permit the addition of an external random force while still obtaining a 
closed-form kinetic equation (Markovian or otherwise) are under investigation, 
but no results have as yet been achieved. 
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VIII. DISCUSSION 


The main point of the article ia considered to be an indication 

that it is possible to derive a systematic kinetic theory of fluid 

turbulence from the Liouville equation for the Fourier coefficients of 

13 

the fluid variables, a line of inquiry begun by Edwards. A kinetic 
equation has been derived and is shown to possess a number of the 
requirements that any reasonable tunetic equation must have: Conservation 

laws, positive-definite spectral densities, and an H-theorom. Some 
value of the kinetic equation obtained may derive from its use in a 
phenomenological description with V^ k relaxation times obtained by 
fitting decay data. 

The major lack in the theory is any reliable information about the 

relaxation predicted by the complicated linear operator H (Eq. (24)). 

There are enough relaxation times v"J k at our disposal, in the absence of 

quantitative estimates of their values, to match virtually any relaxation 

process with the numerical solutions of Eq. (33) . Other quantitative 

. 40 

theories of relaxing two-dimensional turbulence have shown a gross 

insensitivity to the details of the dynamical description, and it may be 

that a similar insensitivity exists here. 

Comparisons with data and with other analytical theories of decaying 

turbulence await numerical investigation of Eqs. (31), (32), and (33). 

We have devoted little attention to the inclusion of viscosity or stirring 

forces, though some calculations of viscous, forced equilibrium solutions 

41 

to the Liouville equation have been reported by Thompson. 
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